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Abstract 

Methods of Lie group analysis of differential equations are extended 
to weak solutions of (linear and nonlinear) PDEs, where the term "weak 
solution" comprises the following settings: 

(a) Distributional solutions. 

(b) Solutions in generalized function algebras. 

(c) Solutions in the sense of association (corresponding to a number of 
weak or integral solution concepts in classical analysis). 

Factorization properties and infinitesimal criteria are developed that allow 
to treat all three settings simultaneously, thereby unifying and extending 
previous work in this area. 

Key words. Algebras of generalized functions, Lie symmetries of differ- 
ential equations, group analysis, Colombeau algebras. 
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1 Introduction 

Local symmetries for equations with weak type solutions, such as, e.g., conser- 
vation laws, involve different constraints depending on the framework in which 
the equations are analyzed. The aim of this paper is to study symmetry prop- 
erties of differential equations involving singular (in particular: distributional) 
terms through an analysis of symmetries in distribution spaces and generalized 
function algebras, as well as associated (i.e., weak type) symmetries. A main 
ingredient in our analysis will be the determination of infinitesimal criteria for 
these solution concepts. 

Investigations in this direction have been initiated by Methe (|^), Tengstrand 
(Q), Szmydt and Ziemian (131;^, 33 1) and have been systematically pursued 
by Berest and Ibragimov 

(H ITbH ll mi) in the distributional setting. More 
recently, in |l| , an extension of the purely distributional methods applied 

so far has been given that allows to also consider nonlinear equations involv- 
ing singularities. The basic tool allowing for such an extension is Colombeau's 
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theory of algebras of generalized functions. In what follows, on the one hand 
we are going to continue the analysis of and on the other hand we shall 
establish connections between the distributional criteria developed in || and 
the Colombeau-type methods given in pl| ]. Moreover, the framework of gener- 
alized function algebras enables us to study symmetry properties of associated 
solutions (i.e. of weak or integral solutions) by the same methods. 
To begin with, let us fix some notations concerning group analysis of differen- 
tial equations and Colombeau's theory of generalized functions. Our principal 
reference for symmetries of differential equations is , whose terminology we 
shall follow closely. Let Ai be an open subset of Rp+' (in what follows, p will 
be the number of independent variables of a system of differential equations, q 
the number of dependent variables) and G a Lie group acting regularly on A4. 
For X £ W, u eW and (x, u) e M we write 

(i, u)^ g - {x, u) = {^g{x, u),^g{x, u)) 

If Sg(a;,u) = Sg(a;), g e G, then G is called projectable. Elements of the Lie 
algebra of G as well as the corresponding vector fields on M will typically be 
denoted by v. 

Vni X = W,U = W. By identifying a function f -.D.^ X with its graph 
T f — {(a;, f{x)) I X G il} C A" X W, the action of 5 G G onto / is defined (locally) 

by 

<?•/ = ./=(*<, ° (id;t X/)) o (Sg o (id;, x/))-i 

where id;i' is the identity mapping on X] in the projectable case this specializes 
to 

5 • / = *s ° (idA- x/) oSg\ 

Set = UxUiX ...xUn, = MxUiX ...xUn- Here Uk = , with 

coordinates Uj, a = I, . . . ,q, J = {ji, . . . ,jk) G {!,... ,p}'' such that ji < j2 < 
■ ■ ■ ^ jk andpfe — (^^fc^^)' ^^"^ number of different partial derivatives of order 
of a scalar valued smooth function of p variables. Elements of W^"-* are denoted 
by u^"\ A^(") is called the n-jet space of M and we set N := dim(7W^"^). The 
coordinates on M^"'^ will also be written as (zi, . . . ,Zp, Zp+i, . . . , z^) = {x, m'")). 
For any / : C A" — > K', the n-th prolongation or n-jet of / is the function 
pr^"^ / : — > formed by / and its derivatives up to order n. The n-th 
prolongation of a group action g or vector field v is written as pi^^^g or pr'"^, 
respectively. 

Let E be a system of differential equations with p variables and q unknown 
functions of the form 

A,(a:,u(")) -0, i = l,...,s (1) 

(with Ai smooth for all i). We shall henceforth assume that ([|) is nondegenerate 
(i.e. locally solvable and of maximal rank, see j2^). Any f : Q C X which 
solves the system on its domain will be called a solution. This amounts to saying 
that the graph of the 71-jet of /, F^,"' is contained in the zero-set Sa of A. A 
symmetry group of S is a local transformation group on X xU such that if / is 
a solution of the system, 5 G G and g ■ f is defined then also g ■ f is a solution 
of E. 
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Next, let us shortly recall some basic definitions from Colombeau's theory of 
generalized functions 0, |jri|, For notational simplicity we 

are going to work in the so-called "special" Colombeau algebra Gi^), defined 
as the quotient algebra £]\i{Q)/N'{Q), where 

£M{n) {(ue)ee/ e C°°(f7)^ : Vif CC 17,Va e 3pe N with 

sup |9"we(a;)| = 0{e-P) as e 0} 

sup |9"Me(a;)| = 0(e«) as e -> 0} 

Here / = (0, 1]. G{^) is a differential algebra (with componentwise operations) 
and H. G{^) is a fine sheaf on M". The equivalence class of (ue)eg/ in G{^) 
will be denoted by J7 = c\[{u^)^^i] or c\[u^] for short. 
We shall make use of the function spaces 

V{n) ^ {/ e C°^{n) : supp(/) compact} 

S{R") = {/ eC°°(M") : > Va e sup(l + < oo} 

xeR" 

Om{^) - {/ e C°^\fl) : Va G Njj' 3p > sup(l + |:E|)-P|9"/(a;)| < oo} 

2?(ri) is the space of test functions on fi, elements of iS(M") and Om(^") are 

called rapidly decreasing and slowly increasing, respectively. 

The algebra Gt{^) — £t{^)/J^t{^) of tempered generalized functions is defined 

by 

£r{n) = e (OM(f^))' : Va G N^' 3p > 

sup(l + |x|)"f = 0{e-P) (e ^ 0)} 

Ur{n) = {K)eG/ e (OA/(f^))' : Va £ N^' 3p > V g > 
sup(l + = 0(e«) (e ^ 0)} 

By componentwise insertion, elements of G{^) and Gri^) can be composed with 
slowly increasing functions. Next, choose p G iS(R") such that / p{x) dx = 1 and 
/ p{x)x°' dx = Q for all a G Ng with |a| > 1. Then £'(r2) (the space of compactly 
supported distributions) is linearly embedded into G{^) via l : u ^ c\[{u*pe)sei] 
(where p^ = prp( j-))- Moreover t coincides with the identical embedding a : f ^ 
cl[(/)eg/] on X'(r2), so V(yi) becomes a subalgebra of G{^) via l. Finally, there 
is a unique sheaf morphism t extending t to C°°( . ) "—^ T^'{-) ^ G{-) (where 
2?' denotes the space of distributions), t commutes with partial derivatives, and 
its restriction to C°° is a sheaf morphism of algebras. The map l defined above 
also provides a linear embedding of iS'(M") into GtO^") commuting with partial 
derivatives and making 

Oc(R") = {/ e C°°(R") : 3p > Va e sup (1 + |a;|)-P|a"/(x)| < oo} 

a faithful subalgebra. So far, Colombeau algebras are the only known differential 
algebras enjoying these optimal embedding properties. Moreover, an intrinsic 
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global formulation of Colombeau's construction on differentiable manifolds re- 
taining all the characteristics of the local theory has recently been achieved 
(p^, p^). Of the various variants of Colombeau algebras we shall also use the 
subalgebra Goo{^) ofQ{Q) consisting of those elements of 5(0) possessing a rep- 
resentative such that supggj |jue||2^oo(Q) < oo. Finally, we shall employ the 
"mixed type" -algebra tJ(]R x fi) whose elements satisfy ^-bounds with respect 
to t G M and tJT--bounds with respect to x ^ fl. 

The ring of constants in Q is denoted by C, its elements are called generalized 
numbers. Elements of C may be used to model infinitesimal numbers (e.g., 
{s)eei is a representative of an infinitely small yet nonzero generalized number), 
which may be viewed as a "nonstandard" aspect of the theory. The support 
of a generalized function U G G{^), supp^ [/, is defined as the complement of 
the largest open subset 57' of fl such that L/|o' — 0. This notion is coherent 
with the embedding 6, i.e. for any T e 'D'{il) we have suppT = suppg(t(r)). 
Finally, we mention the concept of association in the algebra Q: Two elements 
F, G are associated F « G if there exist representatives and of F and G, 
respectively, such that 



lim / {Ge{x) - F,{x))ij{x)dx = 0, for aU V' £ 
Clearly, this definition does not depend on the choice of representatives. The 



concept of association (resp. strong association, cf. 3.2 below) plays a central 
role in the Colombeau framework as in many cases it allows for a distributional 
interpretation of results achieved in G- Particularly in applications to physics 
and numerics it is of fundamental importance (cf. e.g. 0, Q, JSOf). 
With this terminology at hand we can now formulate the main goals of this 
article: Let G be a transformation group acting on the space of independent and 
dependent variables of a system (|l]) of (linear or nonlinear) differential equations. 
We are looking for criteria for G to transform weak solutions of (|l]) into weak 
solutions. More precisely, we shall develop conditions under which G transforms 
I>'-solutions to ©'-solutions (in case (|l|) is linear), ^-solutions to C/-solutions, or 
solutions in the sense of association into solutions in the sense of association 
(in which case (|l|) is to be replaced by A(a;, w'")) w 0). We shall see that these 
questions are in fact closely linked and that criteria for one situation are often 
useful in other cases as well. As an application we study the symmetries of the 
quasilinear hyperbolic system ut + A{u)ux = where A is an s x s matrix with 
C^(M'') entries (cf. ||3^, j^). For treating this system distribution type spaces 
are not convenient, while Q provides a quite satisfactory solution concept (cf. 
[ p7| , Q and the literature cited therein). We consider strongly associated 
solutions and calculate symmetry transformations of such solutions. Also, we 
discuss infinitesimal criteria for symmetries of this system. 



2 Factorization properties, symmetries in Q 

Let G be a projectable group action on some open set C X x U. As 
was already pointed out in the introduction, composition in the framework 
of Colombeau generalized functions requires polynomial growth restrictions on 
the smooth functions. Thus we first single out those group actions which can 
be applied to elements of G (cf. [^). An element g & G with g ■ {x,u) = 
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(Sg(x), u)), is called slowly increasing ii u i-^ ^g{x,u) is slowly increas- 
ing, uniformly for x in compact sets; g is called strictly slowly increasing if 

$3 e Om{M). 

li n C X, U ^ Q{fl) and g is slowly increasing, the action of g on U is the 
element of Q{Eg{n)) given by 



(<f3o(idx«,)oS;i)^ (2) 



Also, in order to be able to insert elements of G{fl) into (|^) we will from now on 
suppose that w^"^ i— *■ A(x, u*^"^) is slowly increasing, uniformly for x in compact 
sets. A symmetry group of (|l|) in is a local transformation group acting on 
X xU such that if [/ is a solution of the system m Q, g ^ G and g - U is defined, 
then also g - U \s a, solution of S '\n Q. 

Let G be a slowly increasing symmetry group of some differential equation 

A(x,u(")) = (3) 

and let U e G{^) be a generalized solution to (|^). Then for any representative 
of U there exists some (n^)^ e such that for all x and all e we have 

A(x,pr(")ue(a;)) = 71^(2:) (4) 

Due to the nontrivial right hand side of (^) it is clear that a direct (component- 
wise) application of classical symmetry methods to Colombeau solutions is not 
feasible. A transfer of classical symmetry groups into the t^-setting therefore 
has to rely on properties of symmetry transformations that are better suited to 
the algebraic structure of Q. The key concept serving this purpose (and, as we 
shall see shortly, at the same time applicable to symmetries of V- and other 
weak solutions) is that of factorization: 

Let G = {g,, I 77 £ R} be a (classical) one parameter symmetry group of (0) 

9v ■ (2;, u) ^ (S^(a;, u), $,,(a;, u)) 

Then by ||8|, Prop. 2.10 there exists a smooth map Q -.W (W open in 

M X A4("), X C W) such that 

A(pr(")g,(z)) = Q(r;,z)A(z) (5) 

Throughout this paper, >V will denote an open set as specified above. By |^ , 
Th. 3.4 for any smooth u : C M-^ ^ M'' such that gr/U exists we have 



A(^^(a;,u(a;)),pr(")(g^u)(^,,(a;,u(a;)))) 
= Q(?7, X, pr(")u(a;)) • A(x, pr("'M(x)) . 



(6) 



In particular, 77 ^ is a symmetry group of (Ij) in Q if (cc, it*^"'') Q^v(ri,x, 
pj-(") u(x)) is slowly increasing with respect to 7?"^ uniformly for x varying in 
compact sets ([^, Proposition 3.5). Thus the need for conditions ensuring that 
the Q^y remain well behaved (in the above sense) arises. Theorem 3.8 of pT[ 
shows that for scalar differential equations possessing a stand alone term (i.e. 

= c for some k > p) this is indeed always the case. Our first aim is to 
generalize this result. 

Since (|l|) is nondegenerate it follows that the Jacobian J(A) of A has rank 
s on the zero-set of A. The following result uses a mild strengthening of this 
assumption to derive a factorization property adapted to the polynomial growth 
restrictions necessary for applying nonlinearities to elements of Q. 
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2.1 Theorem Let G = {g,, | rj G (— ?7oi%)} be a slowly increasing symmetry 
group of system ^ and suppose that there exist p < ki < ■ ■ ■ < kg < N such 

that, setting J^^ ^^^(A) := ( , the following conditions are satisfied: 

(i) z — > (z', A(z)) is injective, where z' — (zi, . . . , zjt^, . . . , z^^, ■ ■ ■ , zn). 

(ii) z (det(Jfcj....^fc^ (A)))~"'^(z) is defined globally and is slowly increasing, 
uniformly for (zi, . . . , Zp) varying in compact sets. 

2 

Then there exists a smooth mapping Q : W — s- K*^ which is slowly increasing in 
z G uniformly for zi, . . . ,Zp varying in compact sets such that ^ holds. 

In particular, (Q) holds for any smooth it : f2 C Rp — > R"? such that grjU exists. 

Proof. Without loss of generality we may suppose {ki, . . . , kg} = {N — s + 
1,...,A^}. We set z' — (zi, . . . , z^r^^), z" = {zn^s+i, z^) and define 
A : ^ by 

A(z)^(z',A(z)) 

Since det(J(A)) = det(Jfcj_...^fc^(A)) 7^ 0, A is a diffeomorphism by (i). More- 
over, 

A o A~^ = (yi, . . . , yAr) ^ {yNs+i, . . . ,yN) = V ^ y" 

Now set /r,(z) = A(pr'^")5,,(z)). By |2|], Proposition 3.7 and our general as- 
sumption on A, fn is slowly increasing in z, uniformly for (zi, . . . , Zp) varying in 
compact sets. Since rj g^i \s a, symmetry group of (|^) wc have oK^^{y) = 
if 2/" = 0. Thus 

/,oA-i(y) = (/„oA-i)(y',Ty") 1^=0 

- t,j-Af,oK-^){y',ry")dr (7) 
= [S',Jk,...MoK-'){y',ry")dT) ■ y" 

Inserting z = A^^(y) into (|^) we arrive at 

fn{^)= Jfc„...,fc.(/,oA-i)(z',TA(z))dT^ •A(z)=:Q(r;,z).A(z) (8) 

In order to establish the claimed growth properties of Q, since and A are 
slowly increasing, by the chain rule it suffices to estimate Jki...ks {^^^)- The de- 
terminant of this Jacobian is precisely {det{Jki,...,ks i^)))^^ , so the claim follows 
from (ii). □ 

2.2 Remarks 

(i) By dropping the growth restrictions on A, G (in particular, allowing for 
G to be nonprojectable and merely supposing that (det( Jfc^^...,fe^ (A)))^^ 
exists globally), the same proof as above still provides the explicit form 
(H) of the factorization property of general symmetry groups of (|]) , which 
reads 

A(pr(")5^(z)) = (/g Jki,...,kAfn ° • • • , ^fei-i, tAi(z), Zfe^+i, ^^-^ 

. . . , Zfc^_i, rA5(z), Zfc^+i, . . . , zjv) dr) ■ A(z) 
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(ii) A necessary and sufficient condition for (ii) in 2.1 is given by 

CC 3C > 3r > Vz e M^") : 



inf 



(zi,...,Zp}£K 



Uz) >C((l + |zp+i|)...(l + |zw|))-'- 



(10) 



(iii) An extensive compilation of sufficient conditions for global injectivity of 
smooth maps (A in our case) can be found in ]29[ |. As an example we 
mention a result of Gale and Nikaido ([^, ch. 3) stating that any F : 
ri R" (with n a rectangular region in R") is injective if its Jacobian is 
a P-matrix on all of Q (i.e. all principal minors are positive). Moreover, 
global invertibility results for Sobolev functions are given in M . 



2.3 Corollary Under the assumptions of Theorem \2. \ , G is a symmetry group 

of (1; m g. 

Proof. Let u e g{VLY be a solution of (0). Then by (|), for any representative 
(ite)eg/ of u we have, taking into account the projectability of g,,, 



A(a:,pr(")(g^M,)(a;)) 



(11) 



= Q(?7, S_^(x), pr(")u,(S_^(x))) • A(S.^(a:), pr(")u,(S_^(a:))) 
Here, A(S_^(a;), pr^")^^ (S_^(x))) is negligible by definition and Theorem ^.1 



shows cl[a: — * (3(?7, S_^(x), pr^"'''Ue(S_^(x)))] to be moderate. It follows that 
(A(a;,pr("^(f;^Mg)(a;)))£g7 is negligible, which precisely means that gr,u is a so- 
lution in g. □ 

2.4 Example Suppose that system (|^) satisfies = CiSik^ for some p < 

ki < ■ ■ ■ < kg < N ai id n onzero constants c,; . Then (|l|) satisfies conditions 
(i) and (ii) of Theorem 2T. In fact, dct( Jfej_...^fc^ (A)) is constant and since the 
underlying domain is convex the above assumption on (|l|) is equivalent with 

X{z) = CiZfe. + Fi{z') {l<i<s) 



with Fi smooth and z' as in 2.1 (i). From this, injectivity of z — > (z', A(z)) is 



immediate. Thus by 2.1 any slowly increasing classical symmetry group of (^ 



is a ^/-symmetry group as well. This observation applies e.g. to the system 



Vt 



uu, = 



considered in |£1[, Example 3.6. Also, Theorem 3.8 of |2l) is a special case of 
this setup (for s = 1). 



3 Associated and distributional symmetry 
groups 

In this section we are going to investigate symmetries of system (|^) in the sense 
of association, i.e. we shall be concerned with group actions that transform 
solutions of 



A,(x,m(")) «0 (z = l,...,s) 



(12) 
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into solutions in the sense of association. Such group actions will be called sym- 
metries in the sense of association or, by slight abuse of terminology, associated 
symmetries. Also, we will give a rather general criterion for transferring clas- 
sical symmetry groups of linear systems to the distributional setting. In what 
follows we will only consider projectable symmetry groups. 

3.1 Definition u = (ui, ...,Uq) € Gi^)"^ is a solution to ([l^) (also called an 
associated solution to (|l|) j if u has a representative (wie, ...jUqi;) € ^mI^)"^ such 
that for every ip G 'D{fl) 

J Ai{x, pr'^'^'>u,{x))Lp{x) dx as e ^ 0, l<i<s (13) 

The set of all associated solutions to ^ will be denoted by Aa- The set of all 
u S {Qoo)'^ satisfying ( |7^ j will be termed AB/\- Let A C Aa- A symmetry group 
77 — > 5^ of ^ will be called an A-symmetry group if gjjU G A whenever U G A 
and g^iU is defined. 

In what follows, C^(f7) (fc e Nq), the space of compactly supported C'^-functions 
on fl is equipped with the inductive limit topology of its subspaces C|^(f^) = 
{/ e C''{n) I supp/ C K} {K compact in Q). 

3.2 Definition Let u — {ui,...,Uq) G Aa{^) (resp. u G ABa{^)) and let 

k 1-1 

k G Nq. u is called a k-strongly associated (^-associated) solution to (Up if it 
has a representative (ui^, ...,Uq^) G EMi^Y such that for every set B C C^{Q) 
which is bounded in C^(il) we have 

lim sup / Aj(a::,pr'"^Ue(a:;))<y9(a::)(i.T = {I < i < s) 

The space of k-strongly associated solutions to (|l|) is denoted by AS a cind we 
set ABSa '■= AS A n Goo- The corresponding symmetry groups are called ASa- 
and ABS\-symmetry groups, respectively. 



3.3 Remark Associated solutions (i.e., solutions of type ^a) play a central 
role in applications to numerics (cf. the remarks in Section |^). Moreover, by 
imposing slight changes in their definition, Colombeau type algebras can be 
adapted to a wide variety of solution types. Thus Goo is particularly useful for 



the investigation of shock wave solutions of conservation laws (cf. |3.4 |3.6| , |3.7 
below), which fall into this class of generalized functions. It will turn out in 
Theorem 3.5 (i) that the uniform bounddedness of solutions is essential for the 
characterization of ylBiSA-symmetry groups. 



3.4 Example Consider the Riemann problem 

ut + = , . 

u{x,0) ~ Ui + (Ur — Ui)H{x) 

where H denotes the Heaviside function. For ui > Ur the unique weak solution 
to this problem is u{x, t) — ui -\- {ur — ui)H{x — ct), {x, t) S K x [0, 00), where 
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{ui+ Ur)/2. Let e V{R), 6* > 0, / 61 = 1 and set H^{y) = /_^^ 0{x) dx, 
5e{y) = -9{-). We are looking for 1-strongly associated solutions to (|lj) of 
the form Ue{x,t) = ui + (ur — ui)H^{x — ct), with c to be determined. Thus 
let i? be a bounded subset of C^(M x (0, oo)). This means that there exists 
K CCM. X (0, oo) such that supp (p K for all cp ^ B and 

sup{\d°'(p{x,t)\ \(pe B, {x,t) G K, \a\ < 1} < oo 

Letip eB. Then 

[uet(a;,t) + u^{x,t)ui,x{x,t)]'ilj{x,t) dxdt = 



I 

Jr 

j 



x(0,oo) 

[-c{ur - Ul)dxH^{x ~ ct) + 

x(0,oo) 

^dx{{ui + {ur - ui)H^{x - ct))'^)]ip{x,t) dxdt 

[-c{Ur - Ui)H^{x - ct) + l-{{2ui{Ur ~ Ul) 
x(0,oo) ^ 

■ H^{x — ct) + {ur — ui)^H^{x — ci))^)]'0a;(a;, t) dx dt 
^ 1 

[c{Ur - Ul) - - {ul - uf)]^{ct,t) dt 

by dominated convergence, uniformly for ip E B. Thus as above is a 1- 
strongly associated solution if and only if c = (u; + Ur)/2, i.e. if and only if the 
Rankine-Hugoniot jump condition is satisfied. 

The transfer of classical symmetry groups into symmetry groups in the sense of 
association is again governed by factorization properties. Sufficient conditions 
for this transfer are provided by the following result: 

3.5 Theorem Let G — {gri \ rj G {^VOtVo)} be a slowly increasing symmetry 
group of admitting a global factorization of the form (jsj). Then 

(i) If Q depends exclusively on rj, x and u then g^i is an ABSA-symmetry 
group of 

(ii) If Q depends exclusively on x and rj then for any k > g,/ is an ^<S^- 
symmetry group of (jlp. 

Proof, (i) Let </) G and u G ABS Ai^)- It follows from (|) that 

u G Q oc C^ri i^))"^ ■ Moreover, by (^), for any 1 < i < s we have 



A^{x,pr^"\g^Ue){x))(f>{x) dx 

= 5l/ (3ji('7,S_^(x),Ue(S_,,(a;)))Aj(S_,,(a;),pr(")ue(S_^(a;)))0(a;)da; 
j 

= V / Aj{x',pY(^"'>u,{x'))Q,j{r,,x',u,{x'))cj){E„{x'))\detDE,,{x')\ dx' 



(*) 
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If (j) varies in a subset of X'(S^(ri)) which is bounded in Cc then due to the fact 
that the are globally bounded, (*) varies in a Cc-bounded subset of I?(ri). 
Hence the above expression tends to zero, uniformly in (j) and e. 
(ii) With (f) as in (i), 

A, (x, pr(") {grjUe){x))(l){x) dx 

= J2 f A,{x',pY^"^u,{x')) Q,,{Tj,x')cb{E^{x'))\detDE^{x')[ dx' 

(*) 

which tends to zero since (*) is C^-bounded if varies in a C^-bounded subset 

ofP(s^(r!)). □ 

3.6 Example Let ^ be an s X s matrix with C^-entries a*-' = a*-' (u^, u^), i,j 
— 1, s on M*. We consider the quasilinear system 

/\{u,Ux,ut) = ut + A{u)ux (15) 

(for solvability resp. unique solvability of ( p^ ) we refer to |30j and |^) and the 
action of a classical symmetry group G ~ {grj \ rj G (^VOtVo)} 

i = 2i(a;,t) 1 

i = 22(a;,t) y 77 e (-r]o,??o) (16) 
u — ^,^{x,t,u) J 

where E^^E^ e C°"(f^),$^ G (C°°(17 x W)Y , rj G (-770,%)- We are going to 
show that G is an ^5^-symmetry group of (|l^) for each k provided that 

(i) du^ri does not depend on w, i.e. G acts linearly on the dependent variables. 

(ii) does not depend on x. 



Denoting as above the Jacobian of by J2^, and Uf are given by the 

first and second column of the matrix 
where ^li(x, t, u, Ux, Ut) is defined by 



first and second column of the matrix valued function $^ G {C°°{Q, x R'^*))''^^, 



{[<^^x{x,t,u) ^^t{x,t,u)\+^^{x,t)[ux ut]^ {JE^{x,t)) ^ (17) 
Moreover, 

A(u,{t£,'Ut) =: \{x,t,u,Ux,ut) 

= / {^){x,t,u,ux,eut + {l-e){^A{u)u^))d9 X {ut- {-A{u)ux)) 
Jo out 

= Qrjix,t,u)A{u,Ux,Ut), 

Note that Qjj does not depend on Ux, Ut by our assumption on G and the explicit 
form of A. Now 



Ar, ^ ^l{x,t,U,Ux,Ut) ^ + A{^r,{x,t,u))^li{x,t,U,Ux,Ut) 



(18) 
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1 



$7,*] + $r,« [mo; Ut]) 



1 



det JS„ 



771 ■ 



det J2„ 



■=■1 $ 



$1 



Thus 



1 



1 



det JS„ 



■=■2 ■ 



det JS^ 



where Is is the s x s identity matrix. This imphes 



Is — A($ ) 

A^(a;, i, m, u^^Ut) = ■=■2 u^;, Ut) = Q,,A. (19) 



^r]x' — '7]t ' — ^Tyt' — 'TjX 



By assumptions (i) and (ii) it follows that Q,, does not depend on u. Thus our 
claim follows from ^ (ii) . 



3.7 Example Let / S C"'^(K) and let F be a primitive of /. The generalized 
Burgers equation 



Ut + f{u)ux = 

u{x, 0) = U; + {Ur - Ui)H{x) 



has a weak solution of the form 



u{x, t)^ui + {ur - ui)H{x ~ ct), {x, t) eWx (0, T) 



(20) 



(21) 



where c — {F{ur) — F(ui))/{ur — ui). (If /' > and ui > u^, then this solution 
is unique.) The same arguments as in the case of B urg ers' equation in 3.4 imply 
that ( pi|) is the 1-strongly associated solution to ( po| ) if and onl y if c is of the 
given form. In this special case of (|l5|), conditions (i) and (ii) of 3^ are in fact 
necessary for shock wave solutions to be transformed into 1-strongly associated 
solu tions. In fact, if (i) or (ii) is violated then from the explicit calculations in 
3.6 it follows that the integrand in the proof of 3.5 (ii) will contain unbounded 
terms. 

More complex factorizations arise in case / is invertible. Indeed, by a straight- 
forward explicit calculation, the infinitesimal generators 

wi = xtd.j, + Pdt + Y(u){^ - f{u)t)du 

W2 - X^dx + Xtdt + f{u)t)du 

give rise to the group actions 
Gi: S:^x{rj) = i = t{r]) = u = f-\7]x + f{u) - T]f{u)t) 
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with factors 

^ (l-r?t)-V' 

/'(/-i(»72: + /(u)-r;/(u)t)) 
(1 - r;:E)3/'(«) 



Q2 



Theorem 3.5 raises the question of finding criteria for classical symmetry groups 
to display the favorable factorization properties used above. As an important 
case where a general result is available we now turn to systems of linear PDEs. 
With a view to applications in distribution theory (cf . section ^) we also restrict 
our attention to group actions which act linearly on the dependent variables. 

3.8 Theorem Suppose that m A is a linear differential operator: 

N 

M2)= alizi,...,Zp)zk + aoizi,...,Zp) [1 < i < s) (22) 



Furthermore, let 5^ be a one parameter symmetry group of which acts 

linearly in the dependent variables: 

g^x, u) = (Erjix), <^^{x) ■ u + -^^{x)) (23) 

// there exist p < ki < ■ ■ ■ < kg < N such that Jk-^ fc3(^) globally non- 
singular then conditions (i) and (ii) of 2A are satisfied and Q in (^) depends 
exclusively on 77 and x = (zi, . . . , Zp). 

Proof. Using the same notations and conventions as in the proof of Theorem 
2.1 we have A(z) = (zi, . . . , zat-s, A{x) ■ z + ao(z)) where the s x N matrix A{x) 
is of the form {A'{x),A"{x)) with A'{x) a,n s x {N - s) matrix and A"{x) = 
>/fci,...,fc3 (A) invertible. 
Thus' 

A-'(y) = (2/',A"(yi,...,yp)-i • [y" -A'{yi,...,yp)-y' -aoiyi,...,yp)]) 

(24) 



is affine linear in yk for k > p. By ( |23| ) and |g8[, (2.18) we have 

pr(")5^(z) = (S^(a;), $^(a;)zp+i + *r,(a;), . . . , Zk, . . . , zn) (25) 

where 

TV 

^k= b'[{v,x)zi + b'^{rj,x) (26) 

with certain smooth functions bf (Actually, the upper limit N in these sums is 
only attained for terms corresponding to highest order derivatives but for the 
following argument only the general form of ( p6| ) is of interest). Hence both 
and A~^ are (affine) linear in Zk for k > p. It follows that the matrix 
Jki,...,ks{fri ° A~^) is independent of Zfc for k > p. This observation, together 
with (@), finishes the proof. □ 
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3.9 Remark The conclusion of Theorem 3.S remains vaHd for a semiHnear 
system 

N 

^ X! o.l.{zi,...,Zp)zk + ao{zi,...,Zp,Zp+i,...,Zp+g) {1 < i < s) 

k=p+l 

(27) 

provided that ki, . . . ,ks correspond to indices of highest order derivatives of the 
dependent variables. This follows immediately from an inspection of the above 
proof (only the form of Oq in (BJ) changes). 



3.10 Corollary Let G be a symmetry group of the linear (resp. semilinear) 
system (22h (resp. (|2^) j such that the assumptions of Theorem ^ (resp. of 
Remark 3^ ) are satisfied. Then G is an AS /^-symmetry group of (2^) (resp. 



Proof. Immediate from 3.5 (ii) and 3.8 
3.11 Example Let 



□ 



A = Mt + a{x, t)ux + a{){x, t, u) — 0, {x, t G R) 

then setting z = {x, t, u, Ux,ut) and using the notation of |3.8| we have s = 1, kg — 
5, A{z) = (zi, . . . ,Z4, Z5 + a{zi,Z2)zi + ao{zi,Z2, Z3)), A"^) = (yi, . . . , ^4, ^5 - 
a{yi,y2)y4 - ao{yi,y2,y3)), and 

pr(l)5^(z) = (S^(zi,Z2),$r,(zi, 2:2)^3 +^'r;(2l,2;2),Ej=3^i('?'^l'^2)Zj + 

btiv, 21, ^2), Ei=3 bjiv, Zi,Z2)zj + 6^(77, zi, Z2)) 
A(prW5,(z)) = E-=3&^.+^o + «-(E,t3&^.+&o) + ao 

Hence 

d 

Q{y) = ^(^'' ° ^ bl{rj,yi,y2) + a{yi,y2)b2ir],yi,y2). 

Turning now to the distributional setting, we first note that the most general 
group actions applicable to distributions are those which are projectable and 
act linearly on the dependent variables, i.e. those which are of the form (p3|). If 
u e T>'{n)'^ then the action of is defined by 

where Sl^ denotes (componentwise) distributional puUback, i.e., 

iriu),^) = {uiy),ip{f-\y))\detDif-'){y)\)} {u £ V'in'),ipe Vm 

for / : il J7' a diffeomorphism. We will sometimes also write u o f instead of 
f*u. 

3.12 Definition Suppose that ^ is linear and let G be a local transformation 
group acting linearly on the dependent variables. G is called a distributional (or 
T>'-) symmetry group of ^ if it transforms distributional solutions of ^ into 
distributional solutions. 



4 INFINITESIMAL CRITERIA 



14 



a. 12 is the most general definition of distributional symmetry groups. More 
restrictive notions (as introduced e.g. in [^) will be discussed in the following 
section. 

3.13 Theorem Let G be a symmetry group of the linear system (^2|) such 
that the assumptions of Theorem 3.8 are satisfied. Then G is a distributional 
symmetry group of (22). 

Proof. By and (^ , for any smooth w : f7 ^ we have 

A(x,pr(")(g,ti)(a;)) = Q(r/,S_,(a;)) • A(S_,(x),pr(")u(S_,(x))) 

Now suppose that u £ P'(ri)'? is a solution to (|l|). Choose some (p S V{W) 
with J 4>{x) dx — 1 and set ul — Ui * (pe, = {ul, . . . ,ul). Then u\ is smooth 
and converges to Ui in T>' for e — > 0. Since g^ acts linearly on the dependent 
variables we also have that w'^^'^g-qUe) — > pr'-^^ff?/'^) Hence 



A(x,pr(")(5,7/)(a:)) = lim g(,7, S_„(x)) • A(S_,(a:), pr(")7.,(S_,(a;))) 
- Q(77, S_,(x)) • A(S_,(a;), pr(")^i(S_,(x))) = 
in V , which concludes the proof. □ 



4 Infinitesimal criteria 

In this section we develop infinitesimal criteria for finding symmetries applicable 
to all three settings of interest (distributional, weak, Colombeau). The results 
introduced here also establish a direct connection of our approach to symmetries 
of distributional and weak solutions to that given in 1^ (cf. also [|[ |[ D). 
In systems of the form 





Ul 








Uq_ 







, or, for short L{x, D)u ~ F 



(28) 



where 



E 

|J|=0 



smooth, 



[n the order of L) and F — x ^ F{x) G I?'(r2), C W , are examined. The form 
( ^8|) of writing a system of linear PDEs (which, of course, is equivalent to (^2|)) 
provides the advantage of allowing to derive very concise forms of infinitesimal 
criteria for factorization properties, cf. (|3C| ) below. 

In Berest's approach, symmetry groups of (|2^ ) are defined via factorization 
properties: Let G — {g^ | r/ G (^ %,'7o)} be a projectable one parameter group 
acting linearly on the dependent variables. Thus g,^{x, u) — (2(77, x), ^{rj, x, u)) 
is of the form (i = 1, . . . ,p, fc, Z = 1, . . . ,q) 
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G is called a symmetry group of (Eq) if there exists a smooth matrix valued map 
{ri,x) Q{ri,x) such that (see (1.6)): 

L{i, D){g^u)(i) ~ F{i) = Q(77, x){L{x, D)u{x) - Fix)) . (29) 

for all u G 'D'(W')'^ Note that this form corresponds exactly to (^) with Q 
depending on 77 and x exclusively. 

(1.7) gives the following infinitesimal criterion for the validity of (p9|): 



[^D, L]u + L{a{x)u + /3(x)) - ^DF ^ {Lu - F) 

or] ri=0 



(30) 



where 



dr] 



>7=0 



a 



dip 
dr] 



r)=0 



9-0 
9r7 



?j=0 



The equivalence of (^9|) and ( pOD raises the question whether factorization prop- 
erties for general systems (|^) of differential equations can always be character- 
ized by infinitesimal conditions similar to (30). For A smooth, an affirmative 
answer is given by the following result, contained implicitly in |^ (cf. 4.2 (ii) 
below) . Since its method of proof will form the basis for our generalizations to 
the G- resp. I?'-settings we state it explicitly. 

4.1 Proposition Let G = {gri \ rj G (— ?70;'7o)} be a one parameter group with 
infinitesimal generator v acting on M. . Then the following are equivalent 



(i) There exists a smooth mapping Q : W R'* such that 

A(pr(")g,(z)) = Q(77, z) ■ A{z) ((77, z) e W) 

(ii) There exists a smooth mapping Q ^ such that 

pr(")v(A)(z) = g(z)-A(z) (zeTW^")) 



(31) 



(32) 



Proof, (i) ^ (ii): Noting that pr(")v (v the infinitesimal generator of r/ g^^ 
is a vectorfield on A^'") whose flow is precisely pr'^"^(7^, differentiation of ( pl| ) 
with respect to r/ at 77 = gives ( ^ (with Q{z) = ^|^_gQ(??j z)). 
(ii) =^ (i) : dH) yields the following hnear ODE for A o pr^") : 



|^(Aopr(")<?,(z)) 
Aopr(")5^(z)|^=o 



Q(pr(")g^(z)).Aopr(")g^(z) 
A(z) 



(33) 



Let Q{r],z) be a principal matrix solution to (33) (i.e. the i-th column of Q is 
precisely the solution with initial value e^). Then we immediately obtain the 
unique solution to ( p3| ) in the form 

Aopr(")5^(z) = g(77,z).A(z) 

□ 



4.2 Remarks 
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(i) From the explicit formulae given in the above proof it follows that Q de- 
pends exclusively on z G M'^''^ for fc < n (e.g. exclusively on (zi, . . . , Zp) = 
x) if and only if the same is true of Q: indeed it suffices to note that by 
[|||, (2.20) 7r^opr(")g^ = pr^'^^g^ (where tt^ : M^"'> M^^'^ is the natural 
projection). 

(ii) By eq. (2.26), for any nondegenerate system (|^), ( |32|) is equivalent 
to pn^v(Ai) = on the zero-set of A (1 < i < s) which in turn (by p8|] , 
Th. 2.71) is necessary and sufficient for v to generate a one-parameter 
symmetry group of (^. Thus (^ ) is precisely the infinitesimal version of 
the global factorization (||). 

(iii) It follows immediately from the definition of pr^"' g,, that ( pT| ) is equivalent 
with 

A(S^(x, m(x)), pr(")g^(u)(S^(a;, u{x)))) , . 

= g(?7,a;,pr(")M(a;)) • A(a;,pr(")u(a;)) Vm G C°°(f2)' Va; £ ^ ' 

where £7 C Rp runs through all open sets. A similar reformulation is valid 
for (^. 



In the Colombeau setting the general form of (Q) allows for A itself to be a 
generalized function. More precisely, we suppose that A G Qt^M^^'^Y . Thereby, 
the admissible symmetry transformations themselves will become generalized 
functions, so-called (projectable) generalized group actions 5 G t/r(K x MY'^i . 
Thus g is supposed to satisfy t/r-bounds with respect to the group parameter 
77 G R and ^-bounds with respect to (x,u) G M {^^, Def. 4.8). Many of the 
infinitesimal methods of classical group analysis can be recovered in this setting. 
For a detailed analysis we refer to The proof of the analogue to 4.1 in the 



present situation requires the following auxiliary result on solutions of linear 
ODEs in g-. 

4.3 Lemma Let A e ((-»/o, ^/o) x K'')™^ such tfiat for all K CC ^ 



sup / \\Ae{s,x)\\ds = 0(\\og{e)\) 
xeK Jr 

Then for each uq G ^/(Rp)™ the initial value problem 

dtu{t,x) ~ A{t,x)u{t,x) ,„ > 

u{0,x) — uo{x) 

has a unique solution u in ^t((~'7o, ?7o) x M^)™. Setting U the matrix with 
columns the unique solutions with initial conditions Ci (1 < i < m), the solution 
to is given by U{t,x)uo{x). We call U a principal matrix solution to ^3^). 

Proof. We only sketch the argument (for details, cf. |l^, ^). Choosing 
representatives (^e)^ of A and {uoe)e of uq, by the corresponding result in 
the C°°-setting we obtain representatives {ug)^, {Ue)e satisfying the claimed 
properties for each fixed e. Then 



- / \\Aeis,X 

Jo 



Ueit,x)\\ <\\UeiO,x)\\+ / \\A^is, x)\\\\Ue{s, x)\\ ds , t G (-??o, ?7o) 
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Thus Gronwall's inequality and the supposed growth restriction on A imply 
moderateness of the representatives and unique solvability. □ 

Note that the corresponding statement for initial value problems in the sense of 
association (i.e. replacing = by « in (|35|)) is false since unique solvability of linear 
ODEs breaks down in that context. As an easy example take Ue{x) = ^ cos(£^a:) 
(x S K). Then m' !=s but u is not associated to any constant. 
An element u of G{^) is called globally of L°° -log-type (cf. if it possesses a 
representative (ue)e with sup^gj^ l'"e(2;)| = 0{\ log(e)|). After these preparations 
we can state (for the notion of C/-n-completeness, see Def. 4.15): 

4.4 Proposition Let A e (/t-(A^("))'' and let g be a Q-n-complete group action 
on A4 . Consider the statements 

(i) There exists Q G (^^(W))"' with 

A(pr(")g,(z)) =Q(77,z)- A(z) m GriWY (36) 

(a) There exists Q e ^(M'"^)''' 

pr(")v(A)(z) = Q(z) • A(z) in 6;(X'"^)' (37) 



Then (i) implies (ii). IfQopA^^^g^ is an element of Qr{Wy satisfying the 
growth property given in ^.c then (ii) implies (i). If Q and Q are supposed to 
he globally of L°° -log-type then (i) and (ii) are equivalent. 



Proof. Using 4.3, the proof proceeds along the lines of 4.1 



□ 



From the pointvalue characterization of Colombeau generalized functions given 
in |21[ it follows that (pq) is equivalent with 



A(S,(x),pr(")g,(u)(S,(x))) 

= Q(?7, X, pr(")M(x)) • A(x, pr(")u(a;)) in g{VLY Vu G Q{nY m 
(cf. ||l|, Lemma 4.13 and Prop. 4.14). 



(38) 



4.5 Remark If in A e C°°((-77o, 770); Oc(R"))'" and uq G Oc(K"), then 
the solution to (|5|), U{t,x)uo{x) belongs to C°°((-77o, 770); Cc(R"))™- Thus, if 
we suppose A e C'c(K^)" in then we have Q e C'^{{-r]Q,r]Q)]Oc(^^)Y\ 
Q e Oci^^)" and (i) and (ii) in ^ are equivalent 

Turning now to the distributional setting, we first note the following result on 
linear ODEs in V. 



4.6 Lemma Let A e 

problem 



|P)" and let a e P'(IRp) 



Then the initial value 



dfU = A ■ u 
u(0, .) — a 



mC°°(R,r>'(MP)™) 
in V'{RP)"^ 



(39) 



has the unique solution Q ■ a where Q is the smooth principal matrix solution of 
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Proof. That Q • a is a solution follows easily by regularizing the initial data via 
convolution with a standard moUifier and then taking the distributional limit 
of the solutions to the resulting smooth problems. Uniqueness follows from 
uniqueness of the corresponding smooth initial value problem by observing that 
for any solution u e C°°(M, X»'(MP)™) of (|3|) with a = we have dt{Q~'^u) = 0. 
But then u = Q ■ c with c a constant vector which is necessarily 0. □ 

Let us suppose that A is of the form (H) with E V'{W) and al e C°°{W) 
{l<i<s,p+l<k< N). This is the most general form of differential 
operators applicable to elements u of V{W). Moreover, we suppose that the 
group action is of the form (23) (also the most general action applicable 
to distributions). We consider A as an element of V^R^) by embedding oq 
as ao (8) Iw-p into 2?'(M^). For the following result, to simplify notations we 
suppose that = R^' x R"? and that g is defined on all of M . 



4.7 Proposition Under the assumptions formulated before 4^6, the following 
are equivalent 

(i) There exists a smooth mapping Q : R x M*' ^ such that 

(pr(")g^)* A = Q{ri, . ) • A m C°°(M, 2?'(A^("))^) (40) 



(ii) There exists a smooth mapping Q:W such that 

pr("V(A)=Q-A mP'(A^("))" 



(41) 



(iii) There exists a smooth mapping Q : R x M-p — > R" such that 

A(S,(.),pr(")(g,w)) = Q(77, .).A(.,pr(")u(.)) inC^{R,V'{Wy) (42) 
for all u e V'{W). 



Proof, (i) <^ (ii) : Using |4.6| , the proof is again identical to that of (for the 
distributional identity ^((pr'^"^5^)*A) = pr("W(A) used in the argument, see 
e.g., [|4|, Th. 3.7). 

(i) 4^ (iii): With the notations introduced in (|2^), (|25|), (|2^ ) and z — {xi, . . . , Xp, 
Zp+i, . . . , zjv) we have 

N / N \ 

A,(pr(")g,(z)) = E H{v,x)zi + lP,{v.x)\ + 

■j=p+i \i=P+i j 

ao(S(77,x)) ® lAr_p (i = l,...,s) 
It follows that ( ^0[ ) can be written in the form 

N / N \ ^ 

E E a]{^(jl,x))h\{T^,x)\zi+ a){E{7ix))bi{T^,x) + 

i=p+i yj=p+i J j=p+i 

s ( ^ 

a^(S(ry,a;)) Iat-p = ^ (5,„(a;) ^ a'^{x)zk + a'^{x) ®In-p 

m=l \ A:=p+1 
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Hence, introducing suitable smooth functions and distributions Qq {k = p + 
1, . . . , N, i — 1, . . . , s), the proof reduces to establishing the equivalence of 

N 
k=p+l 

= 0, (pm,ipn e P(M) , m = 1,. . .,p, n = p+1, . . . ,N 

and 

N 

{ ^ a'fe(77, a::)zfe(x) + dl{x) ® Iat-^, <y9i(a;i) . . . ipp{xp)) = , 

k=p+l 

frneV{R),m^l,...,p, Zfe eC°°(RP), k^p+l,...,N 

This last assertion naturally splits into a (distributional) invariance property 
of Oq Iat-p resp. Oq and a smooth part, both of which are easily seen to be 
equivalent. □ 



Thus ( |29| ) corresponds precisely to ([42|). In particular, the specific form of the 
infinitesimal criterion for the validity of (^9|) follows from an explicit calculation 
of pr(")v(A) in the notation (due to Berest) introduced at the beginning of this 
section. In fact, we have 



j=l k=l \l=l 



= :'S>k(x,u} 



Then by H, Th. 2.36, pr^v = ELi T^-^ + ELi Eij^o "^"^)5n^ 
where 



i—l 



du 

calculated as follows 



(wj j = -a-^). Thus the r-th component of pr(")v(A) = pr^"^\-{L{x, D)u — F) is 



pr(")v E -pr(")v(F^) 

yfc=i|,7i=o / 

q n p q n p 

= E E E^^^^(«:^Vj+E E *fcV/-£6a.F. 

fc=l|J|=Oi=l fc=l |J|=0 i=l 

q n p q n p 

= E E E^^^^(«s'=)-'^-E E E«s'=^'^^^.-'= 

A;^l I J|^0 i^l k^l \ J\^0 

q n q q n q 

+E E E«'/^''"''"' + E E E""'"'/^''"' 

fc=i|.7|=o;=i fc=i|j|=oi=i 

q n p 

+E E ajW'-E^*^*^'- 

fe=l|J|=0 i=l 

= ([CD, L]u + L{x, D){au + f3) - ^DF)^ 
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This calculation, combined with 4.7 and (p2[) provides a rigorous proof of the 



equivalence of (E9) and (pQ). Moreover, 4.1 allows to derive infinitesimal criteria 



for factorization properties even for systems that are not necessarily linear and 
thereby to obtain workable criteria for finding symmetries of weak or Colombeau 
solutions of such systems. For example, let us consider the semilinear system 

L{x,D)u^ F{u) (43) 

where we shall suppose F G Om(K'') (to allow for an insertion of Colombeau 



functions, note however that 4.8 below does not use this assumption). Further- 



more, let us assume that the group action 17^ is of the form (23). Then we have 



4.8 Proposition Under the above assumptions, the following are equivalent: 

(i) There exists a smooth mapping Q : W — s- such that for all u G C°°(ri) 
{n G W) and all x e n 

i(x, D)ig^u)ix) - F{g^u{x)) = Q{r^, x){L{x, D)u{x) + F{u{x))) (44) 

(ii) There exists a smooth mapping Q : A^'") — » R** such that for all u G 
C°°{n) {Q C RP) and all x e n 

[^D, L]u{x) + L{x, D){a{x)u(x) + p(x)) ~ J{F)(u(x))- 

{a{x) ■u{x)+ I3{x)) ^Q{x){L{x,D)u{x) ~ F{u(x))) ^ ' 

(with J{F) the Jacobian of F). 



Proof. By |j]and ^ (iii) it suffices to calculate pr(")v(A) for A = L{x, D)u- 
F{u) as above. Noting that pr("'v(i^) = J{F)[u) • (a • w + /3), the result follows 
exactly as in the above calculation. □ 

In it was shown that a certain splitting of (|3^) is advantageous for estab- 
lishing a connection between determining symmetry groups of PDEs and group 
invariance of the solutions themselves, in particular with a view to determining 
group invariant fundamental solutions of linear systems. In our more general 
setup, we first note that (using obvious abbreviations) setting K = — a 
and h = a - -^\n=o, (|o|) is equivalent to (^) as well as to (^). Also, (|4|) is 
equivalent to (^): 



{-[L,K] + hL)u-{K + h)F + LI3^Q (46) 
[L, K] = hL , {K + H)F = L(3 (47) 

d 

{-[LJ<] + hL)u- J{F){u)au + — QF{u) - J{F){u)P ^ (48) 

orj i)=o 

We first note the following immediate consequence of (pO|): 

4.9 Proposition Let rj — > be a one-parameter symmetry group of the form 
of the homogeneous system L{x,D)u = satisfying (p^/ If u G G{^)'^ is 
a solution, resp. associated solution resp. strongly associated solution then so is 
Ku. □ 
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4.10 Definition Let rj g^i be a slowly increasing one-parameter group. We 

say that u £ Q{^)'^ is Q-invariant under g^j, Ki-invariant, or ^-invariant, re- 

spectively, if grjU — u, g^fU ~ u, or g^jU « u for all rj. If g,^ is of the form M^) 
and u e I?'(ri)' then u is called V -invariant under g,^ if g^u = u in V . 

4.11 Proposition Let rj g^i be of the form (jg^. 

(i) Let u G {C°°Y '"esp. u ^ Q'^ resp. u e (I?')'. A necessary and sufficient 
condition for u to be invariant resp. Q-invariant resp. D' -invariant under 
gri is that Ku equals (3 in C°° resp. in Q resp. in D' . 

k k 

(ii) Then u is ^.-invariant, resp. ^.-invariant, if Ku~ (3 ~ 0, resp. Ku~ (3 ~ 

Proof, (i) We first note that in each of the possible settings we have 

d_ 

drj 



^{grju) = -Ku + (3 



which is immediate from the chain rule. Thus the conditions are necessary. 
Conversely, for fixed u set 

fiv, = - 0(77, x)u{x) - tp{r], x) 

Then f'{ri,x) = S*Q! • f{ri,x) and f{Q,x) — 0. Thus the claim foll ows from 
unique solvability of linear ODEs in each of the respective settings (cf. |4.3| with 
A smooth (hence automatically satisfying the necessary growth restrictions) and 



16D 



(ii) Setting ri.{ri, . ) — :z*{Kui. — (3) we have 

f'^{r],x) = E*a- fe{ri,x) -\-r^{i],x) 
feiO,x) = 

k 

Also, r^{ri, . ) « resp. ~ for all 77 follows from our assumption (by substituting 
for S^(a:) in the respective integrals). With U (77, x) the principal matrix solution 
to the corresponding homogeneous system we obtain the solution to this initial 
value problem in the form 

fe{v,x) ^ U{r],x) / U{r]',xy^r^{r],x)dr]' 
Jo 

Thus for (/3 G P, resp. f B C , B bounded in C^, we have 

fe{v,x)f{x)dx ^ J r^{i],x) (^U{ri,x) J U{i/ ,xy^ dri'(p{x)j -> 
as e ^ 0, resp. this limit is uniform for (p G B. It follows that /e(?7, .) ~ 0, resp. 

fe k 

fs(T], .) ft! I.e., that g,^Ue ~ u^, resp. gr,Ue ~ u^. □ 
Note that the converse assertion in (ii) does not follow since gr/U — u « resp. 

k 

« for all 77 does not imply any information on ^|q.9?7'^ — —Ku^ + [3. 
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4.12 Example We derive infinitesimal criteria for symmetries ( |16| ) of the quasi- 
linear system (llS) whose infinitesimal generators we write as 



Applying 



d7] 



(,{x, t)dx + t{x, t)dt + ip{x, t, u)du 
to JS,, ^ [^rjx ^rit] + ^7ju [ux Ut] , which is just the 



short form of ([TtI), we obtain 



dri 



d 



t;=0 



+ [Ux Ut] — _J^r, = [i'x V't] [Ux Ut] 



r)=0 



Thus, 



d_ 
dr] 
d_ 
drj 



$1 

V 



))=0 



$1 

V 



Ipx + -ipuUx - ^xUx - TxUt 
t/jt + TpuUt - ^tUx - TtUt 

A„: 



Hence by (|1^) we get the following expression for ^ 

(A($^)) [u^ Ut] 



d_ 

dr] 



ri=Q 



ri=0 



m + miiutYi - ^t[ux]i - rt[ut]i + Y^i;^[^ji[u^\i 



d_ 



ri=0 



[A{u)]l {[^xYl + [^uUUxYl - UUxYl - T,[ut]l) = (^t + A{U)^, 
ii'u -Tt- T^A{u)) Ut + 



(^Aiu)i,,, - ^th - ^xAiu) + J2 ^'1^^ 



(where [a]^ is a reminder that a is b x c matrix). Note that in this expression 



the coefficient of ut equals ^ 



))=0 



Thus the determining system {-^ A)^ 
projectable symmetry reads 



»7=0 



drj 



i(=0 



Qn ■ A for an infinitesimal 



il^t + A{u)tj}x 
dA 

A{u)i}u - £,xA{u) - £_tls +22'^' 



du'- 





(tAu - Tils - T.j:A{u)) A{u) 



or 



ipt + ^(w)V'x = 

[A{u),^^]+Y^rg-^ 



i^th - TtA{u)) + (^J, - T^A{u)) A{u) 



4.13 Example We continue to analyze the last determining system in the case 
of a strictly hyperbolic conservation law 

ut + f{u,v), = 0, (49) 
vt+g{u,v)^ = 0. (50) 
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corresponding to (|15|) with 

^ _ /" fv 

[gu gv_ 

with the characteristic values Ai(u,w) < A2(it, w), G iJ, where iJ is open 

in (cf. [|, 0, |l8); for generahzed solutions, see [|6|). 

We denote by = (^i), 4 = (^l) S {C°°{H)Y^^ the characteristic vectors: 
Ar, = A»r„ ijA^xdj, z = 1,2.' 

We will calculate the coefficients of an infinitesimal symmetry v = ^{x,t)dx + 
T{x,t)dt + 4'{x,t,u,v)du + i'{x,t,u,v)dv for (Q), (pO|). By 4.12 we have 



= 



[A,B] + Mu+Mv = (^l2 - TA)t + (ei2 - rA),^ 
where i3=['^" 

Differentiating ( ^ ) with respect to u and v we have 



(51) 
(52) 



"l" 


+ A 


4^XU 







"0" 


+ A 


4^ XV 


1 


'4^xv 



(53) 



{Bt + (t>xAu + i^xA^) 

[Bt + (t>xAu + i>xAy) 

Putting both together we obtain the matrix equation 

Bt + ABx + (jjxAu + tpxAy = 

Taking the derivative of ( |5^ with respect to x and subtracting from ( |53| ) we 
arrive at 

Bt + BxA + i^h - TA)xt + (^2 - tA)xxA = 
or in shorter form 

Mt + M^A = where M = B + {£,xh - t^A) (54) 

Thus we obtain the system (^l|), ( p4| ) which , while not equivalent to (|5T|), 
( p2[ ) (due to the differentiations used in deriving it), considerably facilitates the 
determination of infinitesimal symmetries. 
By Section 16, we have 



(j){x, t, u, v) 
ipix, t, u, v) 



ai {x — tXi {u, v) , u, v)ri{u, v) , ai € C°° (R x H) 



M{x, t, u, v) — Pi{x — t\i{u, v), u, v)tj {u, v) 



(3i e (C°°(M X H)y''\ {x,t,u,v) eR' X H 



In order to determine ^ and t as well as at, (3i — 

M = B + iJ2 - TxA: 

22 22 



i — 1,2, we use 



'^xfv 
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2 2 2 2 

i—l i—1 i—1 i—1 

By resp. solutions of this system in C°° resp. Q determine infinitesimal 
symmetries in each of the respective settings. 

Note that the above considerations also hold for systems of order s > 2 even in 
non-conservative form with the additional assumption that the coefficients of A 
do not depend on x and t. 
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